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Abstract 

It has recently been realized that brane-antibrane annihilation (a possible explanation for ending 
inflation) may result in defect formation, due to the dynamics of the tachyon field. Studies of this 
possibility have generally ignored the interaction of the brane fields with fields in the bulk; recently 
it has been argued that interactions with bulk fields suppress or even eliminate defect formation. 

To investigate the impact of bulk fields on brane defect formation, we construct a toy model 
that captures the essential features of the tachyon condensation with bulk fields. We study the 
structure of defects in this toy model, and simulate their formation and evolution on the lattice. 
We find that, while the energetics and interactions of defects are influenced by the size of the extra 
dimension and the bulk-brane coupling, the bulk-brane coupling does not prevent the formation of 
a defect network. 
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I. INTRODUCTION 



The early universe is the most likely source for experimental input into physics at very 
high energies. Immediately after the Big Bang, energies as high as the Grand Unified (GUT) 
scale may have been explored. It is hoped that signatures of the processes that took place 
at that time will be visible in the universe today, for instance, in primordial perturbations 
or other relics of the very early universe. 

Recently models for the early universe have become less ad-hoc and more inspired from 
String Theory 0. The brane world model is particularly attractive, since it contains the 
ingredients for a successful inflationary universe. The separation of a brane and an antibrane 
can act as an inflaton field, and the interaction of the pair allows one to calculate the 
corresponding inflaton potential. The collision and annihilation of the brane-antibrane pair 
provides a natural mechanism for ending inflation and reheating the universe. 

One important prediction of a large class of brane inflation models is the formation of 
topological defects, typically networks of cosmic strings — it has been argued [j| that brane- 
anti-brane annihilation in the very early universe leads to the formation of cosmic strings 
only [see also [H]. The formation of topological defects in brane-anti-brane annihilation has 
been well studied in String Theory and Field Theory 0, H, 0] , and represents a particularly 
testable potential signature for these models, since the defect network can have quite different 
properties from that of conventional cosmic string networks 

The discovery of explicit mechanisms for the stabilization of the compact extra dimensions 
makes this class of models even more appealing for phenomenology 0. These types of 
models have been mainly used to describe the inflationary universe, and the predictions these 
models make for the cosmological parameters can already be tested against the experimental 
data. At the same time, the mechanism that stabilizes the compactification poses additional 
challenges since it can make the topological defects unstable. |5( 

The formation of cosmic strings in brane-anti-brane annihilation has also been studied 
in the 4D low-energy effective theory, by analyzing the dynamics of the fields living in 
the world- volume of the brane jfj. One usually considers only the complex tachyon field 
or can include the U {!) gauge field that couples to the tachyon. Bulk fields are usually 
neglected. However, it has been argued in [l[ that the inclusion of the bulk fields in the 
process of defect formation substantially changes the dynamics of defect formation, leading 
to a strong suppression of the defect formation rate. This is because the topological defects 
are themselves branes, so they couple to fields described by closed string modes, namely 
the dilaton and the Ramond-Ramond fields which propagate in the entire bulk space-time 
OH 11- 12l|. The energy of the gradients of the RR field along the compactified extra 



dimensions is large, so the formation of defects is strongly suppressed. More precisely, it 
was argued that deSitter quantum fluctuations cannot account for the formation of lower 
dimensional branes since the typical energy of the quantum fluctuations is of order the 
Hubble constant H while the gradients of the fields spreading in the extra dimensions of 
size R are of the order 1/R. During brane inflation H <C 1/R, so the formation of lower 
dimensional branes is suppressed. 

In this paper we try to analyze the mechanism of defect formation when bulk fields are 
involved. We first construct a toy model which captures the essential features of the full 
10-dimensional model that are responsible for the formation of the topological defects. Using 
this toy model we study the classical evolution of the brane and bulk fields, so the deSitter 
quantum fluctuations are responsible only for seeding perturbations which grow during the 
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subsequent evolution. The final configuration of the bulk field gradients is a result of this 
classical evolution. The energy of the system comes from tension of the brane-anti-brane 
pair, which can be modeled as the potential energy of the tachyon field sitting at the top of 
the potential. In our toy model this corresponds to the potential energy of a complex scalar 
living in the worldvolume of the brane-anti-brane pair. 

We will present the SUGRA Lagrangian for the fields which live both inside the world- 
volume of the brane-anti-brane system and the bulk, and we will construct a toy model with 
the same physical features. We will evolve the equations of motion of this toy model on a 
lattice and observe the formation of the defects. We will analyze the effect the coupling to 
the bulk fields has on the formation process as well as the characteristics of the defects once 
they are formed. 

In section 2 we present the relevant supergravity Lagrangian and the toy model we use. 
In section 3 we analyze the equations of motion of the toy model and discuss the ener- 
getics for a defect-anti-defect pair. In section 4 we describe the results we obtained from 
lattice simulations. Finally the conclusions are presented in section 5. Details of the lattice 
implementation of our toy model are left to a technical appendix. 



II. TACHYON CONDENSATION WITH BULK FIELDS 

The usual approach to defect formation during tachyon condensation is to take into 
account only the fields that live in the worldvolume of the decaying non-BPS brane or 
brane-anti-brane pair. This approach seems to be motivated by the fact that the resulting 
lower-dimensional branes formed during the decay are localized inside the worldvolume of 
the parent brane. However, the final state defects are themselves D-branes, and therefore 
couple to bulk RR-fields. One should include the effects of these fields in the defect formation 
process. 

When studying the formation of topological defects by the brane worldvolume fields only, 
one usually considers a single non-BPS brane with the action: 

S = -T p J d p+l x e~+ V(T)y/det (P [G ab + B ab ] + 2vra' [F ab + d a Td b T}) , (2.1) 

where P [G ab + B ab ] represents the pull-back of the bulk metric and NS-NS two-form field on 
the brane and F ab is the field strength of the Abelian world- volume gauge field [HQE3. 



This action was used extensively to study the evolution of the tachyon field in various special 
settings. The most common one is the spatially uniform field in a Friedmann-Robertson- 
Walker universe. In this case one usually sets the NS-NS two-form field and the brane gauge 
field to vanish everywhere, and study the time evolution of the tachyon field and scale factor 
of the universe. The uniform field does not lead to the formation of defects, it behaves like 
a pressureless fluid known as "tachyon matter" [lif . 

When studying the formation of topological defects during tachyon condensation, one 
usually chooses a flat metric, sets the NS-NS field and the brane gauge field to vanish 
everywhere, but chooses a tachyon field that is both time and space dependent. The equation 
of motion for the tachyon field, 

VI + d a Td«T 2(l + d a Td«Tf 2 V(T) y/1 + d a Td a T 
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is non-linear, and one usually solves the equation around a point where T = 0. There 
the space and time dependence can be approximated by a linear space profile with a time- 
dependent slope, T (t, x) ~ u (t) x, and the resulting equation for u (t) can be solved. The 
solution becomes singular in finite time, the occurrence of the singularity marking the forma- 
tion of the topological defect. This result confirms the String Theory calculation in which a 
linear tachyon profile T (x) = ux reproduces the correct tension of a codimension one brane 
in the limit u — > oo. 

The brane gauge field is usually included in the form of a uniform background electric 
field [17] or a constant gauge potential (Wilson line) The case where both the brane 
gauge field and a worldvolume scalar field other than the tachyon are included was studied 
in [l9|. 

This action is highly non-linear and in order to perform a lattice regularization we prefer 
an expression in which the square-root is expanded to quadratic order in the field strengths. 
Also for a single non-BPS brane the tachyon is a real scalar which cannot be minimally 
coupled to the world-volume gauge field. 

It is therefore more convenient to consider the action for a brane-anti-brane pair. In the 
case of a D9 — D9 pair the expanded action involving the complex tachyon coupled to the 
gauge fields living inside each brane is given in Ref. [2J|: 



S = 2T m I d l0 xe-*e- 2 ™' TT 



1 + 8W In (2) D^TD^T + 



(2W) 2 , . 2 (2W) 2 . - 2 (3a' 2 ( y 



(2.3) 



The two gauge fields live in the worldvolume of each brane and the tachyon field couples 
only with one linear combination: 



DpT = dpT-(A+-A-)T. 



(2.4) 



The brane fields couple with the Ramond-Ramond (RR) bulk fields through the Chern- 
Simons coupling also given in Ref. [Hj: 



qDD 
^RR 



T D9 / C A Str e 



(2.5) 



One can expand the exponential above, and the leading order coupling between the brane 
fields and the bulk RR field involves the same linear combination of brane gauge fields that 
couples to the tachyon: 



S 



DD 
RR 



2na'T m / C p -i A (F+ - F~) 



(2.6) 



We see that the orthogonal linear combination, + A~ does not couple to any other fields, 
so we will drop it from the action. Including the kinetic terms for the RR field, the dilaton 
and the metric, the 10-dimensional action describing the decay of the brane-anti-brane pair 
is: 



s = 4 / ^ dl ° x Y 2 * (R + 2 m 2) " Wy F ' 



-2T, 



D9 



d 10 xe-*e 



-2-Kot'TT 



1 + 8W In (2) D»TD,T+ ^ + / } ^ (Pj-F") : 



, (2-7) 
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where F p is the corresponding field strength for the potential C p _i, F p = dC p -\. The C p _i 
will be the only field we consider here, we will not include the dilaton and the metric in the 
toy model we consider. 

III. A TOY MODEL 

We want to create the simplest toy model which captures the important features of the 
formation of defects with bulk fields included, and at the same time is amenable to a lattice 
regularization. This will allow us to follow the evolution to see whether defects form, and 
how the interactions with the bulk fields affect that formation. We want to include the 
minimal field content that will allow us to study the formation of the defects, so we will not 
include the metric and dilaton fields present in the full 10- dimensional model. 

Since we want to study the formation of codimension 2 defects the brane must have 
at least 2 spatial dimensions, so we will choose a 2+1 dimensional brane. The defects that 
form are 0+1-dimensional and the corresponding bulk field that couples to their worldvolume 
must carry a single index. Therefore we will have a vector field living in the bulk. This field 
corresponds to the C p _i RR field in the 10-dimensional model. We choose the bulk to have 
the minimal space dimensionality, 1 space dimension more than the brane. Inside the brane 
we put the same field content as in the full 10-dimensional model, an Abelian gauge field 
corresponding to the linear combination F+ u — F~ u and a complex scalar field corresponding 
to the complex tachyon T. 

Therefore the model consists on a complex scalar charged under a gauge field living 
in a 2+1 dimensional "brane" (see Table |1J). As in the full String Theory case, the fields have 
a Chern-Simons coupling to the 3+1 dimensional bulk gauge field C M . For our purposes it 
is sufficient to keep only the coupling of the brane gauge field to the bulk field and ignore 
the coupling of the scalar, which is higher order in the string coupling, as we discuss in 
Appendix |XJ 





10D model 


toy model 




dimension 


field content 


dimension 


field content 


Brane 


P+l 


A^ — , T 


2+1 


A ii, (j) 


Bulk 


9+1 


Cp-i 


3+1 





TABLE I: Correspondence between the field content and dimensionalities for the two models. 

We choose the complex scalar to be an Abelian Higgs field instead of a tachyon, since 
the tachyon field develops singularities in finite time and the simulation can only run until 
the first defect forms. The Higgs captures the essential features of the tachyon, but is well 
behaved after the defects have formed, so we can follow the evolution. The action of the 
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model is: 1 



c 



I dx 2 dt 
'm 3 



4g 



2 

brane 



D^D^* - V (0) 



/ FAC+ / dx 3 dt 



Mz 



Ma 



^bulk 



H 



We denote by F = dA the field strength of the field A and by H 



(3.1) 

dC the field strength 



of the field C. The scalar covariant derivative is D, 



One can easily check that 



the action is invariant under gauge transformations for each one of the fields A and C. We 
take the standard symmetry breaking form for the scalar potential, 



V(<P) = X{<p*<p-v 2 /2)'' 



(3.2) 



It is straightforward to derive the equations of motion which follow from this action. The 
simplest one to derive is the equation of motion for the scalar field since the scalar <fi does 
not couple directly with the bulk field: 



d^(p - 2iA"d^ - A^cf) - hjWAp + 



dV 
d<p* 



0. 



(3.3) 



One usually chooses the Lorentz gauge d^A^ = to further simplify the equation above. This 
can be convenient when solving for a static, isolated defect, but is not the most convenient 
one when performing lattice simulations. The equations of motion for the gauge fields are 
the usual ones with an extra term coming from the Chern-Simons interaction: 



ua/3 



0. 



- c cs g 2 h 



c a/3v 



ulk" 



F aP 8{z) = 0. 



(3.4) 
(3.5) 



A. Vortex defects 

At vanishing Chern-Simons interaction, c cs = 0, these equations support a topological 
defect, the Nielsen-Olesen vortex j2]J. In any nonsingular gauge, choosing cylindrical co- 
ordinates (p, 9, z) centered on the vortex, the phase of the scalar field <fi winds by 2im, n 
an integer, in going around the vortex. At large radii, the scalar rests in its broken phase 
minimum |0| 2 = f 2 /2, and the gradient energy which the phase change would give rise to is 
exactly canceled by a gauge field, which in a particularly convenient gauge takes the value, 

yfiMr) = ve me , A e (r) = - . (3.6) 

P 

Close to the core of the defect, the gauge field takes a smaller value and the scalar field 
leaves the vacuum manifold, with = at the exact vortex center. The defect core costs 



Here and throughout we use the [ — H — h] metric convention and geometrical normalization for gauge 
fields, so covariant derivatives are of form = — The scalar field charge is chosen to be 1. The 
mass dimensions of the parameters and fields appearing in the Lagrangian are, [A] = [C] = 1, [<fi] = 1/2, 
[ffUl = °' bbraJ = W = 1, and [4J = 0. 
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FIG. 1: Cartoon of how the bulk H electric field, sourced by the F flux at the defect core, spreads 
out across the extra dimension. The dashed line is the brane, the thick bar is the F field strength 
in the defect core, the lines with arrows are the H = VC° electric flux lines, and the dotted line is 
the periodic boundary in the bulk direction. 

finite energy and carries a net magnetic flux of 27m. Far from the vortex, the fields are 
locally in vacuum and all excitations are massive; so the mutual interactions of two vortices 
or a vortex-antivortex pair are exponentially suppressed with separation. 

Examining Eq. ()3.4|) shows why these conclusions will change somewhat at nonzero values 
for the Chern-Simons coupling c cs . Working in the same gauge, we see that only F p g is 
nonzero, so only C° is sourced by the Chern-Simons term. Its equation of motion is, 

V 2 C° - c cs g 2 hulk F 12 5 (z) = . (3.7) 

The component F\2 = F p ^ acts as a surface charge density sourcing the bulk field. 

Our goal is to try to understand the structure and energetics of an isolated defect, or a 
widely separated defect anti-defect pair, when this Chern-Simons interaction is taken into 
account. As a byproduct, we will determine the force between defects, which will no longer 
fall off exponentially. The force will be attractive between defects of opposite winding 
number, but repulsive between defects of the same winding number; our toy model does 
not include the dilaton and the graviton fields, which can provide compensating attractive 
interaction, so we cannot obtain the BPS states that correspond to two identical, parallel, 
branes. 

One can perform a Kaluza-Klein reduction of the 3+1 dimensional bulk theory down to 
2+1 dimensions. This will yield a theory containing a light <fi and A field and a massless C 
field, together with a Kaluza-Klein tower of heavier C field modes. It is possible to solve 
for the defect structure using this effective theory. Since the brane fields are directly the 
fields appearing in the 2+1D effective theory, their equations of motion change the least. 
The equation of motion for is unchanged, and the expression for A^, Eq. (|3.4|) . is changed 
only by the replacing H with a sum over KK states of H from each KK state. The change 
for the bulk field C is slightly more complicated. The KK decomposition is a Fourier series 
expansion, so 8(z) in z space becomes the constant l/2irR in each KK mode. Therefore, the 
equation of motion for the C° field, Eq. ()3.7|h becomes, 

(V 2 - Ml) Cl = C ~^fF 12 , M k \ =Q\ (3.8) 
with m the index of the KK field (the KK wavenumber). 
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We will not attempt to solve these equations in detail near the defect core, but will 
content ourselves with an understanding of the defect appearance at radii larger than the 
symmetry breaking scale, p 3> l/gbrane^, and larger than the radius of the extra dimension, 
p~^> R. In this case, the appearance of M^ k in the equation for leads to an exponential 
decay of C° away from the defect core for all KK states but the zero mode, which instead 
falls according to the 2D Coulomb's law: 

C°(P) = °-^^Hp/Po), (3.9) 
$ = J dOpdpF 12 = jpdOAe, (3.10) 

where $ represents the total magnetic flux associated with the brane defect, and po is a 
constant of integration. Physically, the charge distribution for the C field is a disk lying on 
the brane at the defect core, but the H electric field spreads out from this disk and becomes 
uniform across the extra dimension at a distance of a few R, as illustrated by the cartoon 
in figure [TJ 

Unlike the case of vanishing c cs , the magnetic flux $ will not equal 2ir. There are two ways 
to see this; by considering the equations of motion, or by considering energetics. Looking first 
at the field equations, the requirement that the scalar field phase change by 2im around the 
defect is topological and is not modified by the Chern-Simons term. However, the presence 
of the Chern-Simons term in Eq. (J3.4)) means that either the field strength or the Higgs 
field current must now be nonzero. The defect core is defined as the region where the gauge 
field strength is non-negligible; outside the core we expect F 12 = 0, so only the Higgs field 
gradient term can balance the Chern-Simons term. This will happen if the gauge field is 
smaller than Ag = n/p (its value in the Nielsen-Olesen solution), corresponding to a flux 
which is smaller than $ = 27m, because in this case the derivative and gauge field parts 
of the covariant derivative Dg will not cancel when acting on <\> = e m9 v/\/2. In detail, the 
equation of motion becomes, 

iSbraae (</>^W " <f jW) = C cs9 l ane d p C° , (3.11) 

^{n-A e / P ) = cJ^I- l ~ (3.12) 

P iTXti Z7T p 

where we used Eq. f!3.9j) in the second line. Solving for $, remembering that $ = 2tt Agp, 
we find, 

f( n -*\ = (<M^) • (3.13) 



2tt J V ^ R J 27r 
$ = 2nn 



-i i ^sfbulk 



2tcRv< 



(3.14) 



B. Defect Energetics 

Alternatively, we can understand this result by considering energetics. If the gauge 
field flux does not equal 2tm, there will be a nonvanishing Higgs field gradient energy. 



S 



However, the larger the flux, the larger the C electric field. The vortex is the minimal 
energy configuration subject to the topolocial condition on the phase of 0. Therefore, the 
gauge field will take the value which minimizes the sum of these two energy densities; 



d , 2nR 



, . . D e <t>*D 6 <j> + ^f-d p C°d p C° ) = . (3.15) 



In the first term, Dg<fi = (in/ p — iAg)<p, and <p*<p = v 2 /2 (since otherwise there would be 
an extensive potential energy). In the second term, C° is determined from Eq. (|3.9j) . Also, 
Ag = Q/2irp, as previously. The minimization therefore requires, 



d 



2p 2 \ 2nJ 2gl^ k \(2n)WpJ 



$ = 2irn 



-2 ^2 



i _|_ C cs5 f bulk 

2nRv 2 



-i 



the same answer we obtained directly from the equations of motion. 

This calculation also allows us to find the energy density far from the brane core, which 

is, 

e=\D 0( p\ 2 + ^-\VCr, (3.17) 

^bulk 

which using the above result for $ equals, 

e = ^^, a = ^%. (3.18) 
2p 2 1 + a 2-kRv 2 V ; 

Note several features of this energy density. First and most important, it falls off as 1/p 2 
as one increases the distance p from the core of the defect. Second, as the Chern-Simons 
coupling c cs is changed from to oo, the energy density far from the core changes smoothly 
from 0, the value for a "local" (Abelian Higgs) topological defect, to v 2 n 2 /2p 2 , the value 
for a "global" topological defect (one for a complex, gauge singlet scalar field). The total 
energy of a defect, J d6pdpe, is logarithmically divergent; 

/v 2 7% 2 a a 
d6 pdp —r- — — = nv 2 n 2 — — In p max / p min . (3. 19) 

2p l 1 + a 1 + a 

Here p m \ n ~ l/gbrane^ is a constant of integration which also absorbs the finite contribution 
from the defect core, and p max is the distance where interactions with the fields of other 
defects become important. For instance, for a defect-antidefect pair of separation r, the 
energy is E = 2irv 2 n 2 (a/ [1+a]) lnr/p min . This behavior shows that a defect antidefect pair 
will feel an attractive 1/r interaction, which is proportional to c 2 s when the Chern-Simons 
coupling is small but has a finite limit when it is large. The Chern-Simons term has replaced 
the exponentially weak long-range interaction with one more typical of Coulomb interactions 
in two dimensions. 

One feature of Eq. (J3.19)) is that the energy associated with a defect is never larger than 
the energy of a global defect, that is, a defect without either the brane field A or the bulk 
field C. We have shown this for the large distance contribution, but it is also true in the 
core; for while the bulk field C costs energy, it only does so if the gauge field A takes on a 
nonvanishing flux. The energetics are that the flux is present because it reduces the scalar 
field gradients, and it will always take a value which reduces the scalar gradients more than 
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it increases field strength energies. Therefore the energetics of defect formation are at least 
as favorable as in a global defect model. Also note that the larger the extra bulk dimension, 
the smaller the parameter a and therefore the energy, as the C° field strength spreads out 
over a larger dimension and therefore carries a smaller field strength squared. 

We can now understand why a network of defects is energetically permitted to form in 
the process of brane-antibrane annihilation, which is assumed to start with the tachyon 
field near zero. The analogous situation in the toy model is that the scalar field starts at 
zero, though with small quantum fluctuations away from zero. Consider an isolated defect, 
hypothetically formed after brane-anti-brane annihilation. There is a finite contribution to 
its energy from the defect core region, and a contribution arising from large p, which grows 
like log p. However, inside the same radius, the potential energy initially available from 
the the scalar field grows like the area, p 2 , and this always wins over the logp dependence. 
We can therefore conclude that for a large enough radius, p, there is always enough energy 
available for creating the defect. The reason is that the energy of the final defect is coming 
from the potential energy of the brane scalar field, or the tachyon in the String Theory case, 
and not from the quantum deSitter fluctuations. 

This study has clarified the behavior of defects in the toy model and argued that their 
formation is not obstructed by energetic considerations. We now will verify these conclusions, 
and study the dynamics of the formation and evolution of the system of vortices which forms 
in brane-antibrane annihilation in this toy model, by studying it via nonperturbative lattice 
techniques. 

IV. LATTICE SIMULATION 

The previous section established a toy model and made preliminary investigations of 
the defect structure and energetics in this model. What we really want is to study the 
dynamics of this model and its defect formation in a nonperturbative way. Therefore we 
have implemented the toy model on the lattice. Excluding the Chern-Simons coupling, which 
will be presented in an appendix, the lattice regularization of the brane and bulk action is 
the usual one for an Abelian Higgs model in 2 + 1 dimensions and a Abelian gauge field in 
3 + 1 dimensions. The fact that the fields are Abelian allows us to use the non-compact 
formulation of the lattice action. 

The lattice equations of motion are most easily evolved in temporal gauge, A = and 
C° = 0. This choice is immaterial so long as we concentrate on gauge invariant observ- 
ables. It makes it slightly harder to verify the field configuration for the vortex discussed 
in the previous section, since the two configurations are related by a time-dependent gauge 
transformation. However, we can still compare lattice studies to our results there using only 
gauge invariant quantities such as field strengths. In particular, defects are easily identified 
by the peak of magnetic field strength at their core, along with the peak in energy density. 

Besides studying the static vortex solutions, the lattice implementation allows us to study 
the dynamical evolution of the fields, including the formation and annihilation of the vortices 
and the flow of energy from the brane into the bulk. 

The goal of the simulations is to see whether defects form, and how the defect network 
subsequently evolves if they do form, when the fields are prepared in a state designed to 
mimic the initial state of the tachyon field right after a brane collision. Namely, we begin with 
vanishing gauge fields A, B, d t A, d t B, and with the scalar field at the unstable symmetric 
point = 0. To simulate initial vacuum fluctuations, we add small fluctuations in the 
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FIG. 2: (Color online) Energy in bulk and brane fields as the bulk size is varied. Each simulation 
begins with a well separated brane-antibrane pair. The brane energy is the solid (red) curve 
and the bulk energy is the dotted (green) curve. The abrupt discontinuity is a brane-antibrane 
annihilation event, which occurs later and later as the extra dimension is made larger, since the 
interaction becomes weaker. 

scalar field. These will grow unstably, triggering the phase transition to the state where 
<p*(j) ~ t> 2 /2 almost everywhere. We look for defects as localized bundles of F 12 field strength 
and \D(f>\ 2 gradient energy, and explore their subsequent evolution. Our lattice simulations 
are generally performed using a lattice spacing av 2 = 0.5, with couplings A = gl iaxie = v<1 

and Sbulk = 1- 

We find that defects do indeed form. At vanishing Chern-Simons number, the defects' 
mutual attraction is exponentially weak, and they persist almost indefinitely in a large lattice 
volume. (We use square, periodic lattices. Because of the periodicity of the lattice, the net 
winding number of all defects, Yl n i is always zero, that is, there are always equal numbers 
of defects and anti-defects, so every defect always has a partner which it can annihilate off 
against.) At finite Chern-Simons number, there is long range attraction and the defects are 
observed to annihilate off much faster, though the larger the lattice volume, the longer it 
takes. 

In order to confirm the two effects discussed earlier, we perform simulations in which we 
change each of two "control" parameters, the size of the extra dimension and the strength of 
the Chern-Simons coupling. This way we can see the effect of the size of the extra dimension 
on the energy of the bulk field and the effect of the strength of the Chern-Simons coupling 
on the magnetic flux through the core of the vortex. 

The effect of -Rtmik is illustrated in Figure |2l which shows the evolution of the energy 
density on the brane and in the bulk for a configuration which was initially evolved with 
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FIG. 3: (color online) The magnetic field strength distribution of a defect-antidefect pair. The 
figure shows the field strength distribution for a series of defect-antidefect pairs in theories with 
progressively smaller Chern-Simons interactions. For the leftmost pair, the Chern-Simons term 
almost eliminates the field strength; for the rightmost pair, the Chern-Simons coupling has been 
turned off and the full field strength is obtained. 

strong Hubble damping until there is a single defect-antidefect pair, and then allowed to 
evolve freely. The total energy of the defect pair grows smaller as the extra dimension 
becomes larger. This is expected since a smaller extra dimension imposes larger gradients of 
the bulk field. We also observe that the defects annihilate more slowly if the extra dimension 
is larger, consistent with our analysis of the energetics. (The larger extra dimension lets the 
bulk field spread out more, weakening the Coulomb attraction of the pair.) For the largest 
extra dimensions we explore (64 and 128 across), the defect pair start close enough together 
that the Coulombic interaction is effectively 3D rather than 2D. 

As the two defects approach each other there is an increase in the energy of the brane 
fields and a decrease of the one in the bulk field. The decrease in the bulk field energy is 
due to the fact that the pair of charges at small separation will generate a dipole field at 
large distances and this field has a lower total energy than two widely separated charges. 
The increase in the brane field energy can be thought of as the energy of motion of the two 
defects plus radiated brane field energy as they interact and annihilate. 

Changing the strength of the Chern-Simons coupling tells another story, see Figure 01 
Looking at the defects themselves, we can see the back-reaction of the bulk field on the 
brane fields. Increasing the strength of the coupling results in a reduction of the magnetic 
flux through the core of the defect, in agreement with the considerations of the previous 
section. 

We also look at the energy transfer rate between the brane and the bulk. One of the 
complications of a symmetry breaking transition such as the one we are studying is, that the 
extra energy density associated with the initial field potential energy has to "go somewhere" 
during the simulation. Frequently people carry out such simulations using dissipative dy- 
namics, meaning the extra energy is eaten up by friction. This is justified in our case if the 
Hubble expansion rate is large. Otherwise, for conservative evolution, the energy density 
typically finds its way into the plentiful short wavelength excitations, where it eventually 
becomes equipartitioned. Because there is a brane-bulk interaction in our simulations, the 
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FIG. 4: (color online) Energy transfer between the brane (solid red line) and the bulk (dotted green 
line) for different values of the Chern-Simons coupling. For large or small coupling the transfer 
is slow, so the bulk does not act as an efficient energy sink. A slow flow of energy can allow the 
brane scalar field to homogenize, reducing the final density of defects. 

bulk acts as an energy sink which can be much more efficient than Hubble expansion. If 
the energy of the brane is drained quickly, the fields will have less time to homogenize, 
increasing the density of topological defects. The question of energy leakage is also relevant 
when studying reheating. 

When studying this effect we use a small lattice 32 x 32 x 8 which allows us to evolve 
the system for a much longer time. As expected (see Figure HJ), for a weak coupling the 
energy is transferred slowly to the bulk, so the brane fields have more time available to 
homogenize. However as the strength of the coupling is increased the energy transfer rate 
reaches a maximum, and then it decreases again. We can understand this effect if we take 
into account the fact that the magnetic flux of the brane gauge field acts as the electric 
charge for the bulk gauge field. We saw before that when increasing the strength of the 
coupling the charge increases. However, when the charge becomes too large, it becomes 
energetically expensive for brane magnetic fields to appear at all; therefore a very large 
Chern-Simons interaction simply prevents brane magnetic fields from taking large values, 
turning off the mechanism available for exchange of energy between brane and bulk. 
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V. DISCUSSION AND CONCLUSIONS 



Previous studies of defect formation during tachyon condensation have neglected the role 
played by bulk fields. It was argued jH that when bulk fields are included, the formation 
of topological defects is suppressed or prevented. The argument was that the energy of the 
deSitter fluctuations is of the order the Hubble constant during inflation, while the gradients 
of the fields along the compact dimensions have corresponding energies 1/R ^> H ; so there 
is not enough energy available for defect formation. 

We have constructed a toy model that captures the essential features of the tachyon 
condensation process when both the brane worldvolume gauge field and the bulk RR field 
are included. In this model the gradients of the bulk field are formed as a result of classical 
field evolution and the energy is provided by the tension of the brane-anti-brane pair. Since 
this energy can also be described as the energy of the tachyon field placed at the top of the 
potential, in the toy model this corresponds to the potential energy of an Abelian Higgs field. 
We argued that inside a large enough area around a defect there is always enough energy in 
the form of potential energy of the scalar field such that there is no energetic constraint on 
the formation of topological defects. 

In terms of energetics, the presence of the bulk with a Chern-Simons coupling increases 
the energy cost of a defect. The increase is greatest if the Chern-Simons coupling is large 
and if the extra dimension is small. However, the energy associated with a defect never 
exceeds the value in the global defect model, that is, the value it would take if there were 
no gauge field or if g^rane = 0- The bulk field also changes the long range interactions 
between defects, so that their mutual attraction or repulsion shows power law, rather than 
exponential, dependence on separation. This speeds up the evolution of the defect network, 
so that defect-antidefect pairs find each other and annihilate more efficiently. Again, the 
defect network evolution is never faster than in the global defect model, which is known to 
support a defect network after a symmetry breaking transition. 

Another effect we observe in the lattice simulations of the toy model is the decrease of 
the magnetic flux unit through the vortex core as a function of the Chern-Simons coupling 
constant. This flux also acts as the charge density for the bulk field and if the Chern-Simons 
coupling becomes too large this charge density becomes very small, reducing the energy 
leakage from the brane to the bulk. 

One drawback of the toy model comes from using an Abelian Higgs field instead of a 
tachyon on the brane. We did this so that we could follow the evolution of the defects after 
they have formed while the tachyon develops singularities. At the same time the brane fields 
do not drop out of the dynamics once the defects have formed. 
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APPENDIX A: LEADING ORDER CHERN-SIMONS COUPLING 



Here we expand the Chern-Simons coupling of the brane and bulk fields to second order 
in a'. The full expression is 

SEE = T D9 JCA Str e 2ma ' r , (Al) 

where 

9-p 

C = E CT^..^ A • • • A dx»> . (A2) 

p=odd ^ >' 

The supertrace of the matrix is defined as, 

Str M = Tr Q ^ ) M ' ( A3 ) 
and is the curvature of the superconnection, given by 

ir =(' F+ - TT DT -). (A4) 
\ DT iF~ — TT J K J 

Here we are interested in the particular case of codimension 2 defects, so the relevant 
coupling will be with the C p _i RR-field that couples with the defects. We want to keep 
both the tachyon and the gauge fields non-zero, so we will expand the exponential inside 
the supertrace in powers of a': 

V -1 7 V ; V — 1 y V DT iF~—TT 
| (2rof Tr / 1 \ / iF + -TT DT \ ( iF+-TT DT 



2 VO-iyV DT iF~-TT J \ DT iF~—TT 

+ ... 

= (2ma') [F + - F~] + (A5) 
(2to') 2 / ' / / ' / iTT(F^ - F~) - -(F~ A /- M -F'AF') +.... 



The terms of the type FAF couple with C p _3 and count only in the formation of codimension 
4 defects. Therefore, in the case of a D9 — D9 pair, the important couplings between the 
brane and the bulk fields are 

T D9 J -iC 8 A |(27ria') [F + - F~] + (2na') 2 [DT A DT - iTT (F + - F~)] } . (A6) 

Since we are interested in the simplest model which involves such a coupling, when we 
construct the toy model we keep only the interaction that corresponds to the leading term 



in a', 



(2W) T D9 Jc 8 A [F + - F~] . (A7) 
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APPENDIX B: HIGHER DIMENSIONAL GENERALIZATIONS 



We can extend the toy model we used to a higher dimensional one which will have 
fewer differences with respect to the full string-theory model. The most straightforward 
generalization would be to consider a 3 + 1-dimensional brane and a 9 + 1 dimensional bulk. 
The field content of the worldvolume theory would still be the Abelian Higgs model, but the 
bulk field will now be a rank-2 antisymmetric tensor field, which has the appropriate rank 
to couple to the 2-dimensional world- volume of a string. The Abelian Higgs model in 3 + 1 
dimensions admits stable string-like defects which in our model will be charged under the 
bulk field. The Lagrangian of the model is: 



/ d 3 xdt 

J Ma 



4<? 



brane 



FAC+ / d 9 xdt 

Mi J Mm 



-H 1 



'bulk 

(Bl) 

where as before we denote by F = dA the field strength of the field A and by H = dC the 
field strength of the field C. The equations of motion now become: 

va/3-y 

d,F^ + i^ ra „e {4>D V 4>* - + CcdLoe — Hafiy = , (B2) 

8,H^ - c cs g 2 hulk ^F a(3 5 (z) = . (B3) 

If we now want to study how the solution for a Nielsen-Olesen string is modified by the 
presence of the bulk field, for the brane fields we make the usual ansatz for a string placed 
along the z-direction: far from the string the scalar and gauge fields are 

V2<i>(r) = ve ine , Mr) = z- , (B4) 

P 

with z a rescaling of the c cs = case which remains to be determined. We observe that for the 
bulk field only the C 03 component is sourced by the Chern-Simons term. The entire analysis 
done for the lowest-dimensional toy model can be applied here with the only difference being 
the expression of the masses for the KK modes, wich will be highly dependent on the details 
of the compactification. For the simplest torus, the expression is: 

6 2 
m 



ML = E^" 



=i 



However, when studying the energetics of the defects, only the zero mode gives a potentially 
log-divergent contribution, and the details of the compactification are relevant only to an 
energy density associated with the core. Also, as in the toy model studied in the main text, 
the (global) defect with the field varying by a 2ir phase around the defect core but both 
the F and H field strengths vanishing gives an upper bound on the defect energetics. As in 
the lowest-dimensional toy model we have: 

C 03 (P) = C -^^Hp/ Po ), (B6) 



$ = 



= J d9pdpF 12 = j pd9A e , (B7) 
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with V the (6-dimensiorial) volume of the compact manifold. This leads to the same result 
for the unit of magnetic flux through the vortex core: 

-i 

(B8) 

and for the energetics of the defect; the string tension involves a log, ln(p max /p min ), with 
Pmax the inter-string separation and p min the string core size; 

T = ™V^lnp max /p min , a^^T^lk. (B9) 
1 + a Vv z 

We could also have considered higher dimensional branes wrapped on cycles of the com- 
pact manifold, but again only the zero modes of the brane fields would participate in the 
formation of the defects. 



$ = 27m 



1 + 



Ccs^bulk 



Vv 2 



APPENDIX C: 10 D GLOBAL DEFECTS IN THE ORIGINAL THEORY 

The easiest setup to study the formation of defects in the original theory is to consider 
the action for a brane-anti-brane pair for only the tachyon field in flat space-time, with all 
the other fields turned off, since the lower- dimensional branes are in fact vortices of the 
complex tachyon field. The action is simply: 

S = 2T D9 J d 10 X e ~ 27va ' TT [1 + 8W In (2) D^TD^T] . (CI) 

One can also consider the action for lower- dimensional brane-anti-brane pairs. The equation 
of motion derived from the acion above is: 

d^T - 2-Ka'TdfjTd^T + — ^ = . (C2) 

As in the case of the real tachyon field, one can approximate the profile of the field with a 
linear one as the vortex will form at the place where T = 0, and solve the resulting equation 
for the slope of the profile. The resulting defect formed in the decay of a Dp — Dp pair is a 
Dp — 2 brane. 

In order to undestand this we have to go back to the calculation of the space-time action 
calculated on linear tachyon profiles and estimate the action in the limit of infinite slope. 
The calculation was done in Ref. [2(| and we reproduce the important points here. In the 
case of a D9 — D9 pair the space-time action for a linear tachyon profile is: 

S (y 1 ) = 2T m I dX 10 e~ 2 ™' TT f[ F (na'y 1 ) (C3) 
J i=i 

where T 1 = ^X 1 /\fa' and y 1 = (u 1 ) 2 . The function F has the expression: 
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and in the large argument limit it takes the form: 

F(x)~^. (C5) 
Calculating the action on the profile y 1 — > oo and y 2 — > oo the authors of Ref. [20] obtain: 

S (y 1 ) = 2T m J dX l0 e< yl ^ 2+y2 ^h [ixa'y 1 ) F (na'y 2 ) 

= 2T D9 J rfX 8 y^yjv^Wv^W -> 4vrVT D9 y dX 8 . (C6) 

The result gives the correct tension for a D7 brane, Td7 = (27rvaj Tbg. Regarding the 
RR charge of the vortex, we can estimate it by using the result Eq. (jA5|) in the expression 
of the coupling between brane fields and the bulk RR fields, 

Srr = T m Jc 8 A e~ 2 ™' TT (W) a dT A dT = ^^~T m J C 8 , (C7) 

which again reproduces the correct result for the D7 brane RR charge. 

These results allow us to obtain an upper limit for the density of defects formed, based 
only on energetic considerations. The brane tension is equal to the mass per unit volume for 

the brane and the result Td p -2 = (27rvo7) Tjj p tells us that we can have at most one defect 

on each patch of area (2n\/a'Y = (2irl s ) 2 where l s is the string length. This is a very large 
density and it shows that constraints other than the energetic ones are more important in 
determining the final density of defects. 

In a realistic model we expect that a very important role will be that of the other fields 
that we have neglected so far, namely the dilaton and the graviton. These two fields have 
universally attractive interactions and their presence allows for the existence of BPS states 
in which there is no interaction between identical, parallel, branes. We expect the presence 
of these fields to also change the evolution of the resulting network of defects, since our 
toy model allows for repulsive interactions between same-charge defects, while no repulsive 
interactions (except for very special situations, see are possible in the full String Theory 
model. 



APPENDIX D: LATTICE IMPLEMENTATION 

The implementation of gauge fields on the lattice was first developed by Wilson ,23| and 
is by now standard. To render the number of degrees of freedom finite, the scalar field <fi is 
taken only to reside at a discrete set of points, a cubic (or square) lattice with spacing a. 
Since we will want a finite-timestep update, it is also convenient to define the theory from 
the very beginning on a spacetime lattice, with temporal lattice spacing a t <C a. The gauge 
fields are a connection (rule for parallel transport) and must be defined on the links (lines 
between nearest neighbor lattice sites). Therefore the lattice variable is 

A£(x) = A-dl~aA», (Dl) 

J X 

with (i the unit vector in the \x direction. In a non-Abelian theory it is necessary to treat the 
exponent of this variable, rather than A L itself, as the natural variable, but in an Abelian 
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theory we are free to consider either or exp(— iA£) as the native variable, and we will 
choose to use A^ (the noncompact formulation). Note that the variable A^(x) is centered 
at x + afx/2 

The terms in the action, Eq. become, 



Mz 



a2<it H H 

x&Mz i 



|e lA ^(b{x+ai) — <t>(x)\ 2 



Mz 



2 v-^ |e lA t(j)(x+a t i) - Mx)\' 
a a t 



a? 



FijF lJ 



a 2 a t 



Mz 



EE 



(D2) 
(D3) 



[A[(x) + A^x+ai) - 4 (x+aj) - Aj(x)] 2 ^ 



and similarly for F 0i F 0i and the H 2 terms. These terms are all very standard in the lattice 
community. For the Chern-Simons term, the magnetic field is the sum of A fields going 
around a square, 



Try 
L a 2 



— A 1 (x) - A 1 (x+aj) - A 3 (x) + A 3 (x+ai) 



(D5) 



which, note, is centered at x + (i+j)a/2. Therefore the C field must be averaged over the 
four sites, x, x+ai, x+aj, and x + a{i+j). Further, since C° lives at the half time-step, it 
must be averaged over C{t) and C(t — a t /2); 



Mz 



F 12 C°^a 2 a t Y F 12 x—Yc° 
^ L 8a t ^ L 



xeMz 



(D6) 



The sums for F 01 C 2 and F^C 1 work similarly. 

Alternatively, one can re-arrange this sum in terms of what A fields couple to each C 
field; 



12 

L 



xdMz 



x£M 



3 



(D7) 



Note that F 12 involves a signed sum, see Eq. (|D5|) . so the "middle" lines in the field strength 
squares above cancel off, the contribution associated with one C° field "link" is, 



F 12 C° 



(D8) 



In the continuum, time evolution requires fixing initial field values <ft, A, C and their time 
derivatives Dq<J), Foi, H^. On the spatiotemporal lattice, we fix the links on two initial time 
slices, t = and t = at- The choice is constrained by Gauss' law, which is the extremization 
of the lattice action with respect to a temporal link A^(x) or C®(x). Extremization of the 
action with respect to the and A variables on the t = a t layer fixes the values on the t = 2a t 
layer once we choose values for the temporal links (choose the time dependent gauge), which 
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is most conveniently done by setting = C® = 0. In the c cs = theory, there is a 1-1 
correspondence between equations of motion and t = 2a t variables; variation of A l L (x,a t ) 
determines A l (x, 2a t ) explicitly. At finite Chern-Simons term, the relations between the 
t — at equations of motion and the t = 2at fields are implicit and must be solved iteratively 
by perturbing in c cs . 2 This works because the Chern-Simons term has one derivative, while 
the term has two; so the iteration converges in powers of a t /a. 

We should also verify that the action used is gauge invariant. Changing the gauge at 
point x means shifting 4>(x) —>■ e~ ld (f)(x) and making a change in the gauge fields which will 
keep the <fr field derivative term unchanged. Examining Eq. (jD2|) . we see that this requires, 

4>{x) -> e- w (f)(x) (D9) 
A[{x) -> A[{x)+e (D10) 
A[(x-ai) -> A[(x-ai)-6 (Dll) 

so each link starting at x is shifted by 9 and each link ending at x is shifted by —9. This 
clearly cancels in the field strength, F 12 ; each corner of the square in Eq. (|D5|) has one 
link entering and one link leaving. This ensures the gauge invariance of the term in the 
action, and the gauge invariance to changes in the A field gauge of the F A C term. What 
about C field gauge changes in the Chern-Simons term? 




FIG. 5: Magnetic fields contributing to the variation of the action with respect to a gauge change 
in the C field at the center of the cube. This sum of A field links exactly cancels-a lattice 
implementation of a "boundary of a boundary." 

To see that the Chern-Simons term is invariant to C field gauge change, first note what 
such a gauge change adds to the action. Each C link entering or leaving the site x is shifted 
by ±9 (+ if it leaves the site). This contributes 9 times a signed sum of F£ v "squares" to 
the action. But when we look at the signed sum in detail, see Figure |5J we find that it is 



2 Specifically, the Chern-Simons term means, for instance, that the variation of the action with re- 
spect to A 1 L (x,a t ) depends on A^(x,2a t ), which is to be determined. One starts by guessing that 
A? (x, 2dt) = (x, at) to evaluate the equation of motion at x, a t . This gives a first guess for A^ (x, 2at) 
(and simultaneously, all other links at time 2at). This guess is used to compute the Chern-Simons term, 
redetermining the equation of motion at time a t and a better guess for the variables at 2a t . The process 
is iterated until numerical convergence, typically of order 4-7 iterations. 
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the surface of a box, such that each link appears in two F^ u } with opposite orientation; the 
sum cancels. 



[1] G. Dvali and A. Vilenkin, "Formation and evolution of cosmic D-strings," JCAP 0403, 010 

(2004) |hep-th/0312007] . 
[2] G. R. Dvali and S. H. H. Tye, "Brane inflation," Phys. Lett. B 450, 72 (1999) 

|hep-ph/9812483| . 

S. H. S. Alexander, "Inflation from D- anti-D brane annihilation," Phys. Rev. D 65, 023507 



(2002) hep-th/0105032 . 

A. Mazumdar, S. Panda and A. Perez-Lorenzana, "Assisted inflation via tachyon condensa- 
tion," Nucl. Phys. B 614, 101 (2001) |hep-ph/ 01 070581 . 

C. P. Burgess, M. Majumdar, D. Nolte, F. Quevedo, G. Rajesh and R. J. Zhang, "The 

inflationary brane-antibrane universe," JHEP 0107, 047 (200 1) |hep-t h/0105204]. 

G. R. Dvali, Q. Shafi and S. Solganik, "D-brane inflation," pIep^t h/01052o'3T~ 

C. Herdeiro, S. Hirano and R. Kallosh, late "String theory and hybrid inflation / acceleration," 

JHEP 0112, 027 (2001) |hep-th/0110"27"T] . 

B. s. Kyae and Q. Shafi, "Branes and inflationary cosmology," Phys. Lett. B 526, 379 (2002) 



hep-ph/0miMj. 



C. P. Burgess, P. Martineau, F. Quevedo, G. Rajesh and R. J. Zhang, "Brane antibrane 
inflation in orbifold and orientifold models," JHEP 0203, 052 (2002) |hep-th/01li"025| . 
J. Garcia-Bellido, R. Rabadan and F. Zamora, "Inflationary scenarios from branes at angles," 
JHEP 0201, 036 (2002) |hep-th/0112"l47] . 

R. Blumenhagen, B. Kors, D. Lust and T. Ott, "Hybrid inflation in intersecting brane worlds," 
Nucl. Phys. B 641, 235 (2002) |hep-th/0202"T24] . 

K. Dasgupta, C. Herdeiro, S. Hirano and R. Kallosh, "D3/D7 inflationary model and M- 
theory," Phys. Rev. D 65, 126002 (2002) |hep-th/0203"019l . 

N. Jones, H. Stoica and S. H. H. Tye, "Brane interaction as the origin of inflation," JHEP 
0207, 051 (2002) |hep-th/0203l"63"l . 

M. Gomez-Reino and I. Zavala, "Recombination of intersecting D-branes and cosmological 
inflation," JHEP 0209, 020 (2002) |hep-th/0207"278l . 

S. Kachru, R. Kallosh, A. Linde and S. P. Trivedi, "De Sitter vacua in string theory," Phys. 
Rev. D 68, 046005 (2003) |hep-th/0301240] ; 

S. Kachru, R. Kallosh, A. Linde, J. Maldacena, L. McAllister and S. P. Trivedi, "Towards 
inflation in string theory," JCAP 0310, 013 (2003) jhepjth/0308055). 

J. P. Hsu, R. Kallosh and S. Prokushkin, "On brane inflation with volume stabilization," 
JCAP 0312, 009 (2003) |hep-th/ 03 11077] . 

H. Firouzjahi and S. H. H. Tye, "Closer towards inflation in string theory," hep-th/03 12020. 
E. Halyo, "Inflation on fractional branes: D-brane inflation as D-term inflation," JHEP 0407, 
080 (2004) lhep^th/0312042"]. 



E. Halyo, "D-brane inflation on conifolds," hep-th/0402155 . 

C. P. Burgess, J. M. Cline, H. Stoica and F. Quevedo, "Inflation in realistic D-brane models," 
|hep-th/0403"TT9l 

A. Buchel and A. Ghodsi, "Braneworld inflation," hep-th/0404151 



[3] S. Sarangi and S. H. H. Tye, "Cosmic string production towards the end of brane inflation," 



21 



Phys. Lett. B 536, 185 (2002) hep-th/0204074 . 

N. T. Jones, H. Stoica and S. H. H. Tye, "The production, spectrum and evolution of cosmic 

strings in brane inflation," Phys. Lett. B 563, 6 (2003) |hep-th /0303269|. 
[4] N. Barnaby, A. Berndsen, J. M. Cline and H. Stoica, "Overproduction of cosmic superstrings," 

JHEP 0506, 075 (2005) |hep-th/ 04 12095] . 
[5] L. Pogosian, S. H. H. Tye, I. Wasserman and M. Wyman, "Observational constraints on cosmic 

string production during brane inflation," Phys. Rev. D 68, 023506 (2003) hep-th/ 0304188 . 

G. Dvali, R. Kallosh and A. Van Proeyen, "D-term strings," JHEP 0401, 035 (2004) 

|hep-th/0312005| . 

E. J. Copeland, R. C. Myers and J. Polchinski, "Cosmic F- and D-strings," JHEP 0406, 013 
(2004) |hep-th/0312067] . 

E. Halyo, "Cosmic D-term strings as wrapped D3 branes," JHEP 0403, 047 (2004) 
|hep-th/0312268| . 

L. Leblond and S. H. H. Tye, "Stability of Dl-strings inside a D3-brane," JHEP 0403, 055 
(2004) |hep-th/0402072] . 

K. Dasgupta, J. P. Hsu, R. Kallosh, A. Linde and M. Zagermann, "D3/D7 brane inflation and 
semilocal strings," JHEP 0408, 030 (2004) |hep-th/0405247l . 
[6] P. Brax, J. Mourad and D. A. Steer, "Tachyon kinks on non BPS D-branes," Phys. Lett. B 
575, 115 (2003) |hep-th/0304197l . 

E. J. Copeland, P. M. Saffin and D. A. Steer, "Singular tachyon kinks from regular profiles," 
Phys. Rev. D 68, 065013 (2003) |hep-th/0306"294"l . 

P. Brax, J. Mourad and D. A. Steer, "On tachyon kinks from the DBI action," hep-th/0310079. 
[7] J. A. Minahan and B. Zwiebach, "Effective tachyon dynamics in superstring theory," JHEP 
0103, 038 (2001) |hep-th/0009246| . 

K. Hashimoto and S. Nagaoka, "Realization of brane descent relations in effective theories," 
Phys. Rev. D 66, 026001 (2002) |hep-th/0202079"l . 
[8] I. R. Klebanov and M. J. Strassler, "Supergravity and a confining gauge theory: Duality 
cascades and xSB-resolution of naked singularities," JHEP 0008, 052 (2000) hep-th/0007191| . 
S. B. Giddings, S. Kachru and J. Polchinski, "Hierarchies from fluxes in string compactifica- 
tions," Phys. Rev. D 66, 106006 (2002) |hep-th/0105097] . 

S. Kachru, R. Kallosh, A. Linde and S. P. Trivedi, "De Sitter vacua in string theory," Phys. 

Rev. D 68, 046005 (2003) [hep-th/0301240l . 
[9] J. Polchinski, "String theory. Vol. 1: An introduction to the bosonic string," 
[10] J. Polchinski, "String theory. Vol. 2: Superstring theory and beyond," 
[11] C. V. Johnson, "D-brane primer," |h"ep-th/0007170l 
[12] M. R. Douglas, "Branes within branes," |hep-th/9512077l 

[13] M. R. Garousi, "Tachyon couplings on nonBPS D-branes and Dirac-Born-Infeld action," Nucl. 
Phys. B 584, 284 (2000) |hep-th/0003"T22| . 

E. A. Bergshoeff, M. de Roo, T. C. de Wit, E. Eyras and S. Panda, "T-duality and actions 
for nonBPS D-branes," JHEP 0005, 009 (2000) |hep-th/0003"22l"] . 

J. Kluson, "Proposal for non-BPS D-brane action," Phys. Rev. D 62, 126003 (2000) 
|hep-th/0004106| . 

A. Sen, "Tachyon matter," JHEP 0207, 065 (2002) ||hep-th/0203265| . 
[14] D. Kutasov and V. Niarchos, "Tachyon effective actions in open string theory," Nucl. Phys. 
B 666, 56 (2003) |hep-th/0304045] . 

M. Smedback, "On effective actions for the bosonic tachyon," JHEP 0311, 067 (2003) 



22 



|hep-th /0310138|. 

[15] A. Sen, "Rolling tachyon," JHEP 0204, 048 (2002) [hep-th/ 0203211 . 

N. Lambert, H. Liu and J. Maldacena, "Closed strings from decaying D-branes," 
|hep-th/0303139[ 

F. Larsen, A. Naqvi and S. Terashima, "Rolling tachyons and decaying branes," JHEP 0302, 
039 (2003) |hep-th/0212248] . 

[16] A. Sen, "Tachyon matter," JHEP 0207, 065 (2002) |hep-th/0203265] . 

A. Sen, "Field theory of tachyon matter," Mod. Phys. Lett. A 17, 1797 (2002) 
|hep-th/0204143| . 

G. W. Gibbons, "Cosmological evolution of the rolling tachyon," Phys. Lett. B 537, 1 (2002) 



hep-th/0204008 . 



G. Shiu and I. Wasserman, "Cosmological constraints on tachyon matter," Phys. Lett. B 541, 
6 (2002) |hep-th/0205003l . 

G. Shiu, S. H. H. Tye and I. Wasserman, "Rolling tachyon in brane world cosmology from 

superstring field theory," Phys. Rev. D 67, 083517 (2003) |hep-th/0207lT9] . 
[17] P. Mukhopadhyay and A. Sen, "Decay of unstable D-branes with electric field," JHEP 0211, 

047 (2002) |hep-th/0208142] . 
[18] K. Hashimoto, "Dynamical decay of brane-antibrane and dielectric brane," JHEP 0207, 035 

(2002) arXiv:hep-th/0204203| . 
[19] T. Mehen and B. Wecht, "Gauge fields and scalars in rolling tachyon backgrounds," JHEP 

0302, 058 (2003) plep^t h/02062l2|. 
[20] P. Kraus and F. Larsen, "Boundary string field theory of the DD-bar system," Phys. Rev. D 

63, 106004 (2001) |hep-th/0012198| . 
[21] H. B. Nielsen and P. Olesen, "Vortex-Line Models For Dual Strings," Nucl. Phys. B 61, 45 

(1973). 

[22] G. Gibbons, K. Hashimoto and P. Yi, "Tachyon condensates, Carrollian contraction of Lorentz 
group, and fundamental strings," JHEP 0209, 061 (2002) |hep-th/0209034] . 
K. Hashimoto and S. Terashima, "Brane decay and death of open strings," JHEP 0406, 048 
(2004) |hep-th/0404237] . 

[23] K. G. Wilson, "Confinement of quarks," Phys. Rev. D 10, 2445 (1974). 



23 



